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Abstract. Using a large sample of pure, slow, short lived K° mesons collected with KLOE detector at 
DA$NE, we have measured the Ks lifetime. From a fit to the proper time distribution we find r(Ks) = 
, (89.562 ± 0.029 s tat ± 0.043 sys t) ps. This is the most precise measurement today in good agreement with 

» I ■ the world average derived from previous measurements. We observe no dependence of the lifetime on the 

' direction of the Ks- 
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1 Introduction 

We have collected very large samples, O(10 9 events), of 
slow A-mesons of well known momentum, with the KLOE 
detector at DA$NE. Kaons originate from the decay of cj>- 
mcsons produced in e + e~ collisions. We have used the 
above samples to measure many properties of kaons such 
as masses, branching ratios and lifetimes, refs. 1 through 
7. The ultimate motivation was the determination of the 
quark mixing parameter V us , see ref. 8. KLOE had not 
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however attempted to measure the Ks lifetime. We pre- 
sent a precise measurement of the Ks lifetime based on a 
sample of about 20 million Ks— >ti + -k~ decays correspond- 
ing to an e + e~ integrated luminosity of 0.4 fb _1 . 

The reaction chain e + e~— ></>, <p—> K l (unobserved) Ks, 
Ks~ >-7T + 7r~, with p0 = 13 McV in the horizontal plane, is 
geometrically and kinematically overdetcrmined. We can 
therefore, event by event, determine the As-meson vector 
momentum px, the kaon production point P and its decay 
point D. From px , P and D we obtain the decay proper 
time of the Ks- A fit to the proper time distribution gives 
the ifs-meson lifetime. The vast available statistics allows 
us to select some 20 million A'g— s-7r + 7r~ decays with favor- 
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able configuration to provide the most accurate and least 
biased measurement of time. Averaging over the sample 
gives a statistical accuracy of ~2 //m in the measurement 
of the kaon mean decay length. For consistency we use our 
value of the kaon mass, Mr:=(497.583 ±0.021) MeV, ref. 
7. 

The KLOE detector has been described in all the ref- 
erences mentioned above, see also refs. 9, 10, 11, 12. In 
particular ref. 8 summarizes the use of the KLOE detec- 
tor in collecting kaon data and reconstructing all decay 
channels. 



2 Data reduction 

Data were collected in 2004 with the KLOE detector at 
DA$NE, the Frascati (^factory. DA$NE is an e + e" col- 
lider operating at a center of mass energy - s /s~1020 MeV, 
the 0-meson mass. Beams collide at an angle of 7T-0.025 
rad. For each run of about 2 hours, we measure the CM 
energy y/s, and the average position of the beams in- 
teraction point P using Bhabha scattering events. Data 
are combined into 34 run periods each corresponding to 
an integrated luminosity of about 15 pb _1 . For each run 
set, we generate a sample of Monte Carlo (MC) events 
of ^3x equivalent statistics. We use a coordinate system 
with the 2-axis along the bisector of the external angle of 
the e + e~ beams, the so called beam axis, the y-axis point- 
ing upwards and the x-axis toward the collider center. 
Ks~ >tt + tt~ decays are reconstructed from two opposite 
sign tracks which must intersect at a point D with rp> < 10 
cm and \zp>\ < 20 cm, where x = y = z= is the e + e~ aver- 
age collision point. The invariant mass of the two tracks, 
assumed to be pions, must satisfy \M„ n — M K o \ <5 MeV. 
D is taken as the decay point. The kaon momentum px 
can be obtained from the sum of the pion momenta and 
also from the kaon direction with respect to the known, 
fixed <j> momentum p^. We call the latter value p' K . The 
magnitude of the two values of the kaon momentum must 
agree to within 10 MeV. If the two tracks intersect in more 
than one point satisfying the above requirements, the one 
closest to the origin is retained as the Ks decay point. We 
refer to the finding of D as vertexing. 

The above procedure selects a Ks— >ir + ir~ sample al- 
most 100% pure. For each event we need the kaon produc- 
tion point P. In fact only the z-coordinate of P is required 
since the interaction region is 2-3 cm long while the other 
dimensions are negligible and the x, y coordinates well 
known. P lies on the beam axis and is taken as the point 
of closest approach to the Ks path as determined by the 
7r + 7r~ tracks. The resolution in zp is about 2 mm. Events 
with \zp \ > 2 cm are rejected. From the length of PD and 
p' K we compute the proper time in units of a reference 
value of t, the lifetime value used in our MC, t=89.53 ps. 
Its distribution is shown in fig. 1 top, histogram a. The 
distribution has an rms spread of 0.86 r and is not sym- 
metric. Time resolution can be improved discarding events 
with poor vertexing resolution. From MC we observe that 
bad vertex reconstruction is correlated with large values 
of Ap = px — Pk, the difference in magnitude of px and 
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Fig. 1. Top. Monte Carlo. Reconstructed time resolution, his- 
togram a, in units of r for the initial Ks sample (rms spread 
~0.86t) and after cuts and geometrical fit, b, (rms spread 
~0.32 r). Bottom. Distribution of Ap—pK —p'pc for data (dots), 
and Monte Carlo (line). The tail at left is due to the initial state 
radiation. 



p' K . Fig. 1 bottom shows the Ap distribution for data and 
MC. We therefore retain events with coso^ < —0.87, 
0.5 < \a^ +K | < 2.2 rad, \M^ - M K \< 2 MeV and events 
with —0.5 < cos0(7r ± ) < +0.5. a W7r is the opening an- 
gle of the pion pair. The definition of a^ +K is slightly 
more complicated. Information about the angle between 
the positive pion and the kaon at the decay point D is re- 
quired. We must also distinguish between the two tt + -k~ 
'V configurations illustrated in fig. 2. Calling r and s the 
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Fig. 2. The two configurations for a Ks—tft Tr decay. 



projections of kaon and positive pion on the {x, y} plane, 
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a 4- „ is defined as 

a„+ K = sign ((r x s) z ) arccos (^—^- 

The angle ct^+ K is defined in {— it, it}. Positive sign cor- 
responds to the configuration of fig. 2, left. All angles are 
in the laboratory system. 

After applying the cuts above, only ~l/3 of the events 
survive while the rms time spread is reduced to 0.63 r. 
Another significant improvement is obtained performing 
a geometrical fit of each event to obtain the production 
point P and the decay point D. We chose a new point 
P' on the beam axis and a new decay point D' on a line 
through P', parallel to the kaon path, so as to minimize 
the x 2 function 



rD' - r D 



(z P i - 



The proper time distribution, after all cuts and the fit, 
is shown in fig. 1 top, curve b. The rms spread in t is 
0.32 r. We check the correctness of the Ks direction using 
a sample of if^-mesons reaching the calorimeter, where 
they are detected by nuclear interactions. The Kl inter- 
action point in the calorimeter together with the known cf> 
momentum gives the Ks direction with good resolution. 
Comparison with the Ks direction as obtained from pi- 
ons shows a negligible difference. The final efficiency for 
Ks^>-tt + tt~ detection is shown in fig. 3 as a function of 
proper time. 
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Fig. 3. Monte Carlo: final efficiency, averaged over all Ks 
directions, as a function of the proper time in a single run 
period. 



The average efficiency depends on the Ks direction, 
is almost flat and in average is ~9%. Errors in the re- 
construction of the pion tracks can bias the position of P 
and D. In fact, the value of Ks lifetime differs by ~6% 
for events with a^+ K > and < 0, where the sign distin- 
guishes the topologies of the di-pion 'V. see fig. 2. 

We do correct for this effect. From MC we obtain the 
correction, A£k, to be applied to the Ks decay length, as 
a function of Ap. The correction is applied event by event 
to the data. The procedure is repeated for each run period. 
After applying this correction the 6% difference mentioned 
above is reduced to ~ 10~ 3 , although the average result is 
only ~2cr (0.1%) different from the result before applying 
it. 



3 Proper time distribution fit 

MC and data, see fig. 1 top, studies show that the time 
resolution is well described by the sum of two Gaussians. 
We write the resolution function, normalized to unity, as 



r(t,r, ai,a 2 ,a) 
a 

exp 



ctiV2tt 
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and the decay function, for a lifetime r, as: 
d(t) = - x exp (--) x 

The expected decay curve, normalized to unity, is given 
by the convolution 



d(rj) r(t — rf) dr/. 



Allowance must be still be made for small mistakes in the 
reconstruction of the decay and production position, D 
and P. A shift S in the proper time is therefore introduced. 
Thus the function which we use for fitting the observed 
distribution is 

f(t, r, oi, <7 2 , a, 5) = g(t - 6). 

The four parameters, a±, 02, a, S in f(t) depend on co- 
latitude and azimuth, and <p, of the kaon and it is not 
realistic to attempt to obtain them from MC. We divide 
the data in a 20x18 grid in cos 9, 4> and fit each data set 
for the lifetime r with the above parameters free. In order 
to improve the result stability, we retain only events with 
I cos#| < 0.5 and < </>< 360°, discarding in this way only 
^8% of the events. We therefore perform 180 independent 
fits only to events in a 10x18 grid. The fit range, —1 to 
6.5 r, is divided in 15 proper time bins. The kaon lifetime 
is obtained as the weighted average of the 180 Tj values 



t(K s ) = {t)=Y, 
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The corresponding x 2 value is x 2 = J2i( T i ~ ( T )) 2 / fj2 ( T i)- 
We find x7 dof = 202/179 for a confidence level, CL, of 
11.4%. The normalized residuals of the 180 fit values Ti 
have an rms spread of 1.1. Tab. 1 gives the average cor- 
relations between fit parameters and fig. 4 top shows a fit 
example. 
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Table 1. Correlation of fit parameters (averaged values). 
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Fig. 4. Top: an example of the fit with x 2 /dof= 8/10, dots are 
data, the line is the fit result. Bottom: proper time resolution, 
in units of t, as a function of <j> and cos#. 



The resolution (y/acrf + (1 - o)o|) versus {0,(f>} is 
shown in fig. 4 bottom. The resolution varies from 0.22 r 
to ^0.27r over the accepted {cosO, 4>} range with an aver- 
age of 0.24 t. The Si values show a dependence on <f> with 
period 2ir corresponding to a shift of the position of P of 
~ — 10 /im in y and ~50 /mi in x. In addition, a very small, 
10~ 4 , eccentricity of the drift chamber is evident. All these 
effects are consistent with mechanical and surveying inac- 
curacies. To ensure that the lifetime evaluation is correct 
to the 10~ 4 level, we correct the value of (r) obtained 
above by the factor f/T^ c =1.00036±0.00019, where t^ c 
is the result of fitting the MC data with the procedure 
described above. 



4 Systematics and result 

Changes in analysis cuts and FV corresponding to a 
~±60% change in efficiency result in a lifetime shift of 
0.024 ps. Varying the fit range gives a shift of 0.012 ps. As 
mentioned in Sec. 2, we use the KLOE value of the kaon 
mass in the kinematic determination of the Ks momen- 
tum and the calculation of [3k- The measurement of fix 



and the decay position are independent. The uncertainty 
on the calibration of p' K gives an uncertainty of 0.033 ps. 
The uncertainty due to Ks mass is 0.004 ps. All fits are 
then performed assuming uniform efficiency versus proper 
time, resulting in an uncertainty of 0.005 ps. Table 2 sum- 
marizes all systematic errors. The result is stable across 



source 


absolute value (ps) 


cuts & FV 


0.024 


fit range 


0.012 


p' K calibration 


0.033 


kaon mass 


0.004 


efficiency 


0.005 


total 


0.043 



Table 2. Systematic error contributions. 



the entire data taking period. As said before, without ap- 
plying the vertex correction the result still remains within 
2 cr of the final result, but stability with the run period is 
lost. Our result for Ks lifetime is: 

t(K s ) = 89.562 ± 0.029 sta t ± 0.043 syst ps. (1) 

Subdividing the data in 9 intervals and summing over 
cos 8 the dependence of the lifetime becomes quite ob- 
vious. The average of the 9 r(ifg) values are of course ex- 
actly as eq. 1 but x 2 /dof=24/8 for a CL of ~0.2%. Enlarg- 
ing the statistical error by a factor y / 24/8 restores x 2 =8 
(CL=43%) and corresponds to r(A' s )=89.562±0.050, an 
error very close to a/0.029 2 + 0.043 2 =0.052 confirming our 
estimate of the systematic error in eq. 1. 

The result of eq. 1 is in agreement with recent measure- 
ments, ref. 13-15, as shown in fig. 5. Including the present 
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Fig. 5. Recent Ks lifetime measurements. The green band 
represents the new world average, r(iCs)=89.567=L0.039 ps. 

measurement, the new world average for the Ks lifetime 
is t s =89.567±0.039 ps, with x 2 /dof = 0.5/3, or CL-92%. 

In KLOE we can measure the lifetime for kaons travel- 
ing in different directions. We choose three orthogonal di- 
rections, the first being {£i, b±} = {264°, 48°} in galactic 
coordinates. This is the direction of the dipole anisotropy 
of the cosmic microwave background (CMB), ref. 16. The 
other two directions are taken as {£2, b 2 } = {174°, 0°} 
and {£3, b 3 } = {264°, -42°}. After transforming the kaon 
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momentum to the above systems, we retain only events 
with inside a cone of 30° opening angle, parallel (+) 
and antiparallel (— ) to the chosen directions and evaluate 
the kaon lifetime. The 6 results are consistent with eq. 1. 
Defining the asymmetry A — (jg — Tg )/ (r^ +Tg), we 
obtain the results of tab. 3. Systematic errors are strongly 



{£, b} 


A x 10 3 


{264, 48} 


-0.2 ±1.0 


{174, 0} 


0.2±1.0 


{264,-42} 


0.0±0.9 



Table 3. Observed asymmetry. Errors are dominated by statis- 
tics. 
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reduced when evaluating the asymmetry. Results in tab. 3 
show all the asymmetries values are well consistent with 
zero. 

A further check has been performed using all KLOE 
data sample (about 2 fb _1 ). The result for the asymmetry 
in the direction of CMB anisotropy, consistent with that 
given m tab. 3, is (-0.13 ± 0.40 sta t) x 10" 3 . No estimate 
of systematic error has been performed. 
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